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We study the effects of an on-site disorder potential in a gas of spinor (spin-1) ultracold atoms 
loaded in an optical lattice corresponding to both ferromagnetic and antiferromagnetic spin de¬ 
pendent interactions. Starting with a disordered spinor Bose-Hubbard model (SBHM) on a two 
dimensional square lattice, we observe the appearance of a Bose glass phase using the fraction of 
the lattice sites having finite superfluid order parameter and non integer local densities as an indi¬ 
cator. A precise distinction between three different types of phases namely, superfiuid (SF), Mott 
insulator (MI) and Bose glass (BG) is done via a percolation analysis thereby demonstrating that 
a reliable enumeration of phases is possible at particular values of the parameters of the SBHM. 
Finally we present the phase diagram based on the above information for both antiferromagnetic 
and ferromagnetic interactions. 


I. INTRODUCTION 

The physical realization of spinor Bose gases is severely 
impeded by largely the instability of the internal atomic 
states. The excited states typically decay through spon¬ 
taneous emission on time scales faster than that what is 
needed for the gas to equilibriate via collisions. However, 
owing to the huge technological success, dilute atomic gases 
with internal degrees of freedom have been realized. Op¬ 
tically trapped spin-1 is among the first ones to be 

reported [1-6] . The progress has been rapid since then 
with the discovery of spin-1 [5, 7] and spin-2 manifold of 
[6], [8-10]. A flurry of activities have taken place 
since. It is beyond the scope of the present work to furnish 
the large volume of literature that exists on the progress 
of the study of spinor Bose gases. There is an extensive 
review by Stamper-Kurn and Ueda on the subject [11]. 

An interesting offshoot of the study of phase diagram 
of spinor Bose gases is to investigate the role of disorder 
present therein. Since inclusion of disorder in the context of 
an atomic gas loaded in an optical lattice is experimentally 
possible via speckle or periodic potentials [12-16], investi¬ 
gation of its influence on the phase diagram is a pertinent 
task. 

Further, with the advent of techniques that can probe 
single atoms, the accessibility of the quantum world has 
been revolutionized, which can not be obtained by statisti¬ 
cal sampling [17, 18]. This has made possible imaging of the 
insulating phase of the interacting atomic species of bosons 
on an optical lattice with a single atom and single site res¬ 
olution. Among other things, such studies will furnish the 
population density of bosons (or the fluctuation therein) in 
different phases and their local tunneling dynamics in an 
optical lattice. 

All of the above discussion would gain importance in the 
context of studying a system of spinor bosons in an inhomo¬ 
geneous environment, that is in presence of disorder. The 
phase diagram of the spinor particles is in general richer 
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owing to the different signs of the spin dependent interac¬ 
tions that render a polar or a ferromagnetic nature for the 
superfluid phase [4], [19, 20]. 

The 'theorem of inclusions’ [21] forbids a direct MI to 
a SF transition in presence of disorder. A BG phase al¬ 
ways intervenes in between. Since a BG phase originates 
due to arbitrarily large, although rarely occurring clusters 
of one type of phase present amidst another, and as sta¬ 
tistically rare events are difficult to sample, elegant meth¬ 
ods such as stochastic mean-field theory [22-24], quantum 
Monte Carlo(QMG) etc [24], [25-27] were formulated and 
used to study the BG phase. 

In this work we employ a single site mean field theory 
(SMFT) [28] which is very extensively and most commonly 
used tool and apply it for studying a Bose-Hubbard model 
(BHM) that recognizes the spin degree of freedom of the 
atomic gas (we call it SBHM). Corresponding to the dis¬ 
order free scenario, the nature of phase transition in the 
antiferromagnetic and ferromagnetic cases are ascertained 
via the average order parameter and the local density fluc¬ 
tuation, that is the compressibility. The nature of the phase 
transition for the odd (occupation density 1, 3, 5..) and the 
even (occupation density 2, 4, 6..) Mott lobes are different 
for the polar fluid, while this is untrue for the ferromag¬ 
netic case ( see discussion on antiferromagnetic and ferro¬ 
magnetic cases in the next section). 

In presence of disorder, we define a simple physical quan¬ 
tity, X defined as the fraction of the sites with finite SF 
order parameter and non integer occupation density that 
will be used as indicators for the MI, BG or SF phases [29]. 

A percolation of the insulating sites is seen to occur with 
patches of the superfluid order surviving, thereby main¬ 
taining a zero macroscopic superfluid order, although com¬ 
pressible. The percolation scenario is robustly tested and 
a percolation threshold, implying the onset of a SF phase 
is obtained with the help of Hoshen-Kopelman(HK) algo¬ 
rithm [30]. Equipped with all of the above, we map out 
a phase diagram for the {F =1) spinor bosons in a two 
dimensional square (optical) lattice, where we specifically 
focus on the emergence and sustenance of the BG phase in 
both antiferromagnetic and ferromagnetic cases. 

We organize our paper as follows. In the next section, we 
briefly review the mean field theory for the SBHM in the in- 
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homogeneous case (in presence of disorder). In section III, 
we present the results for the average SF order parameter 
and the compressibility in both the antiferromagnetic and 
ferromagnetic cases. Further, we study the variation of x 
in presence of disorder and the finite size scaling properties. 
We obtain the phase diagram based on this information and 
present it thereafter. Finally in section IV we outline our 
concluding remarks. 


II. MODEL 


The behavior of ultracold atoms loaded in an optical lat¬ 
tice with hyperfine spin F=1 can be well described by a 
SBHM [31] as, 

H = -t -b h.c) + - e,)ni 

<ij> cr i 


where t is the tunneling matrix, {i,j) are the nearest neigh¬ 
bour sites, /i is the chemical potential, Uo describes the spin 
independent on-site interaction and U 2 is spin dependent 
interaction which arises due to the difference in scattering 
lengths oq and 02 corresponding to S=0 and S=2 channels 
and C/q, U 2 are related to the scattering length by, Uq = 
{Airh?/M)[[ao + 2a2)/3) and U 2 = {A-kTi^/M){{ a 2 - ao)/3). 
The spin dependent interaction, U 2 can have either pos¬ 
itive or negative signs depending upon whether 02 > oq 
or 02 < oo- Thus U 2 is antiferromagnetic (where the SF 
phase is polar) for 02 > oq and ferromagnetic (where the 
SF phase being ferromagnetic) for 02 < oq. The total spin 
at site i is Si = aJ^Fo-cr/Oicr' where Fo-o-' are the components 
of spin-1 matrices and (t=-|- 1, 0, -1. The particle number 
operator, fii = hia = a\^ha where al^(dia) is the 

boson creation (annihilation) operator at a site i. Ci is the 
on-site disorder is introduced by the parameter i, which 
is randomly chosen from a box distribution extended over 
[—A, A] where A is the strength of the disorder and an 
important parameter of our work. The components of the 
spin operator are given by. 


A _ A ^ ^ A A _ A A 

^iz - ^ 2 + ^2 — 

Siy = -^(a|Qai+ — d|_|_aio — SJo®*- 


dl_dio) 


— 2Tli-\-TliQ hi-|- 2Tli(j Ui— hi_j_ 


- 2hi+hi_ -b hi_ -b 2a]’+a]’_aig -b 2di+di-d 


iO 


( 2 ) 


To solve Eq.(l), we shall employ SMFT that decouples the 
system Hamiltonian into sum of the single site Hamiltoni¬ 
ans. The hopping part of the Hamiltonian can be decoupled 
as [19, 32], 

= {dt)dja + at(a]„) - (at)(a]^) -b Saljaj^{3) 


where (O) denotes the equilibrium value of an operator, O 
and 6aja = dja — {dja) denotes the fluctuations in dia-- Now 
introducing the superfluid order parameter at site i as, 

= (^jcr) “ \/^S(T^cr — {dia) — (4) 

where 'ipi=i/^= ^-|- and risa is the su¬ 
perfluid density and ^a is a normalized spinor obeying 

It is obvious that all spinors are related to each other by 
a gauge transformation, e*® and spin rotations, «(«, /3, t) = 
g-iSj,ag-iS„/3g-iS^T (a,/3,T) are the Euler angles. 

The behavior of the SF order parameter is different in the 
antiferromagnetic and ferromagnetic regimes and are enu¬ 
merated in the following [33-35]. 

(i) Antiferromagnetic case: when U 2 > 0, the spinor, ^a is 
given by. 


e = 




' —-^e ''°‘sin(3\ 


003/3 


V2 


e"'" sin/3 


( 5 ) 


If a=P=^ then 'i/j+='ip- ^ 0, ■i/'o = 0 or a=[0, 27r] and /3=0 
or TT then '0+='0_=O, ^ 0. 

(ii) Ferromagnetic case: when U 2 < 0, ^a is given by, 

/l\ / e“*“cos^/3/2 \ 

C 0 = e*® K/2 cosP/2 sm/3/2 (6) 

yoy y e*“sin ^/?/2 y 


li /3 = t: 12 and a = 0 then i/>+ = %/>-, i/jq = '/2ip+. 

Now neglecting the quadratic fluctuations and substitut¬ 
ing the superfluid order parameter in Eq.(l), the SBHM 
can be written as a sum of single site Hamiltonians as. 


< ^ 1 . 


MF 


where 


^ = -t +* E 

a . 

- (^ - ei)ni -b -^niidi - 1) -b -^(S^ - 2n^) 

where 4>ia = Fjj V’icr- Tti® sum j includes all nearest neigh¬ 
bors at the site f of a square lattice with z = A, z being the 
coordination number. 

Our main task is to find the SF order parameter, i/^ia by 
diagonalizing the mean field Hamiltonian in Eq.(7) in 

the site occupation number basis, \nia)- For that, we first 
have to find the matrix elements of , which is done 
via, 

(h.+ ,h.o,n.-|fff^|n'+,n'o,h') = + (8) 

where the hf part includes the matrix elements coming from 
the diagonal part of the mean field Hamiltonian and reads 
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as, 


hi = 4>iai’i'T - ifJ-- ei)ini+ + mo + m-) 

+ ^{ni+ + n^o + m-){ni+ + mo + nt- - 1) 

H- — (271^+71.^0 + 271i+71iO + 'ni+ + 271^0 + + ^i+ 


- 2{m+ + ni+ + 7ii+)) 


,n 7 _j_ ^TiiQ^n'^Q ^rii— ,n^_ 


(9) 


While the off diagonal part is, 


= t 


1 + \/^i+<5ni+,n'^ + l^nio.'n.Jg,n<_ 

“ t “ 4 ^ i 0 ,'^^0 — 1 — 

“1“ ^i—\/hli—^ni+^n'.,^niQ^n'.Q^ni—^n'._—l “t“ 


( 10 ) 


+ 


t/2 


2 1 


\/n^o{nio - l){m+ + l)(nj- + 1) 


,n7 ,+1 h-^ 


Hence we diagonalize the matrix in Eq.( 8 ) to obtain the 
ground state energy Eg{'tpic,) and the eigenfunction 
starting with some some guess value for tpicr- Now from the 
updated wave function 'i>g{4’ia), we compute the SF order 
parameter using, 


V^zcr — ^ g{}pia^\hia\^ 


( 11 ) 


Using this new value of ipia, we again reconstruct the ma¬ 
trix and repeat the diagonalization procedure until self con¬ 
sistency condition is reached [28]. 

The variation of the equilibrium ground state en¬ 
ergy, Eg with occupation number, n shows that it 
almost stabilizes corresponding to n=7, for which 
(ni+,nio, 7 ii-|i 7 j^^| 7 i'+,n'o, 77 .') is a 120 x 120 matrix [see 
Fig.l], The dimension is obtained by considering all pos¬ 
sible combination of 77 ,^+, mo, m- such that, m^ + mo + 
m-=m. The local density, pi and compressibility, m can 
also be computed using. 

Pi = (lE'g(7/’icr)|77.icr|lE'g(7/’icr)); Hi = (Pi) ~ (pi) (12) 


In the clean state, that is when A/Uo=0, the mean field 
Hamiltonian is homogeneous and hence the SF order pa¬ 
rameter is uniform throughout the entire lattice. In the 
strong interaction limit (C/q ^ t), the SF order parame¬ 
ter vanishes indicating that the system is in the MI regime 
where the atoms are forbidden from tunneling into neigh¬ 
bouring sites and thereby are localized to individual sites 
with a fixed number of bosons per site. Thus the MI phase 
is characterized by zero SF order parameter and zero com¬ 
pressibility with a finite gap in the particle hole excitation 
spectrum. 

In the weak interaction limit {t ^ C/q), the SF order 
parameter shows finite value whence the atoms can easily 



FIG. 1. The variation of the equilibrium ground state energy. 
Eg with occupation number n. 


tunnel between the nearest neighbour sites. In the super¬ 
fluid phase, atoms are delocalized over the entire lattice 
and the ground state is a coherent state, which in experi¬ 
ments is seen to form a constructive interference pattern in 
the momentum space. Thus the superfluid phase is char¬ 
acterized by non zero order parameter, non integer boson 
density and finite compressibility with no particle hole ex¬ 
citation gap in the spectrum. 

In presence of disorder, the mean field Hamiltonian be¬ 
comes inhomogeneous and hence the SF order parameter 
will vary from one site to another. Further, in the pure 
case, the system can make direct phase transition from MI 
to SF phase controlled by the system parameter t and Uo- 
But as soon as disorder is present, the BG phase always 
intervenes between the SF and MI phases. The BG phase 
is defined as the zero SF order parameter with no gap in 
the particle hole excitation spectra. Thus the three types 
of phases can be characterized as, 

(i) SF phase: tjji 7 ^ Oj Pi y^integer and > 0 

(ii) BG phase: t/jj = 0, pi =/= integer and m > 0 

(iii) MI phase: tf^i = Oj Pi =integer and Ki = 0 

In order to see the effects of disorder, we first solve our 
mean field Hamiltonian, Eq.(7) self consistently on a two 
dimensional (2D) square lattice of size L x L where we 
have considered T = 30 throughout our calculations (un¬ 
less mentioned otherwise) to obtain the averaged value of 
7/; and p as, 


i/j = 


2 = 1 




- sample 


and K = 


I ^ 
2=1 


A^2 


sample 


(13) 


where sample in the subscript refer to the fact that results 
are averaged over different realizations of disorder. Here 
we have used a maximum of 40 different disorder realiza¬ 
tions and confirmed that the fluctuations are satisfactorily 
equilibriated. 
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III. RESULTS 

In this section we shall present our results by nu¬ 
merically solving the mean field Hamiltonian, Eq.(8) in 
the zt — ^ plane while keeping ?7o=l and considered 
different values of spin dependent interaction U 2 - Since 
U 2 /Uo={a 2 — ao)/(ao -I- 202 ), so U 2 /U 0 extends from -1 
when 02=0 to 0.5 when ao=0. In our work, we choose 
C^2/C^o=0.1 (less than the maximum value 0.5) to study 
the odd and even MI lobes. It may be noted that at 
U 2 /U 0 > 0.5, the odd MI lobes are completely consumed 
by the even MI lobes. For the ferromagnetic case, any value 
oi U 2 in the interval [-1, 0] is seen to yield qualitatively 
similar results and we choose U 2 /Uo=- 0 . 2 . 

The other important parameter is the disorder strength, 
A/C/q. Since the MI phase is characterized by finite 
energy gap in their energy spectrum, in the atomic limit 
(t=0) the width for the MI lobe is energy gap between 
the upper and lower values of the chemical potential, fj, 
corresponding to the particle and hole excitations [36]. In 
the antiferromagntic case, the width of the odd MI lobes 
is Eg=fi+ — fi-=Uo — 2 U 2 , while for the even MI lobes, 
it is Eg=Uo + 2 U 2 [37]. In order to go from a gapped MI 
phase to a gappless BG phase, the disorder strength, A/Uq 
should be greater or equal to the width of the respective 
MI lobes. Since we are using a box disorder from -A to 
A, so the critical disorder strength, Ac should extend 
from -Egl2 to +Eg/2 that is Ac/t/o=0.4 for the odd MI 
lobes and Ac/t/o=0.6 for the even MI lobes corresponding 
to U2/Uo=0.1. Similarly for the ferromagnetic case, the 
energy gap corresponds to Eg=Uo + U 2 and hence the 
critical value is Ac/C/o=0.4 corresponding to C/2/t/o=-0.2. 
Using this parameter values, we shall present our results 
systematically in the following. 

A. The behaviour of the SF order parameter 
and compressibility 

In this section we present our results of the averaged 
order parameter, ip and compressibility, k to characterize 
the MI, BG and SF phases as a function of tunneling 
strength zt/Uo. In the pure case, the ip and R undergo 
a transition from zero to finite values indicating a direct 
MI-SF transition. That is, the system remains in the MI 
phase as long as the tunneling strength, t is well below the 
critical tunneling strength tc^ where the latter is given by 
[38], 

zE = (1/3) (Go + 2t/2)[(2p + 3) - VV + 12p] (14) 

While in the disordered case, the BG phase appears and 
it tries to displace the MI phase making inroads for itself. 
The behaviours of ip and R in both antiferromagnetic and 
ferromagnetic regions is discussed below. 

(i) Antiferromagnetic case: The spin dependent in¬ 
teraction, U 2 for is antiferromagnetic since the 

atomic scattering length 02 > oq. For t/ 2 /t/o > 0, the 
nature of MI-SF phase transition is either of second or 
first order corresponding to the MI lobes with odd or even 
local densities respectively, while for a scalar Bose gas, it 


is always a second order phase transition. 

In the Mott insulating phase, with an odd value for 
the local density per site, the SF order parameter, ip and 
the compressibility, R changes continuously from zero, 
thereby indicating a second order MI-SF phase transition. 
However, the MI phase with each site having even local 
density, the ip and R show finite jumps indicating the first 
order MI-SF phase transition [19], [39]. 

We have also found that the three SF order parameters 
assume the values, such as ipi+=ipi- ^ 0 and ipio=0 in the 
MI phase making the ground state only functions of ipi+ 
and ipi- and hence the equilibrium value of the total spin 
is (S)^=0, which is responsible for making the SF phase a 
polar state [19], [20]. 

In order to visualize the different kinds of phase 
transition corresponding to the MI lobes with different 
occupation densities, we plot ip and R in the pure and 
disordered cases. We choose the values of ^ in such a way 
so that we are at the MI lobes with local density pi=odd 
and even. At /i/[/o=0.4, we are in the vicinity of the odd 
MI lobe with the local density pi = 1 while for p/Uo = 1.4, 
the local density comes out as pi=2 implying we are at the 
even MI lobe. 

In the pure case, we found that corresponding to 
U2/Uo=0.1, the system remains in the Mott insulating 
phase with ip and R as zero till zt/Uo=0.16 for the first odd 
MI lobe at /r/17o=0.4 and zt/Uo=0.26 for the first even 
MI lobe at p/Uo=lA, both of which lie below the critical 
tunneling strength ztc/Uo. Beyond this critical value, the 
system goes to a superfluid phase with non integer local 
densities, p and finite values of ip [see Fig. 2 (solid (black) 
lines)]. 

We shall now include the on-site disorder, and see 
how it affects the averaged SF order parameter and com¬ 
pressibility. The inclusion of can equivalently be treated 
as a site dependent chemical potential, pi {pi=p — Ci). 
Since the SF order parameter and local densities become 
explicitly site dependent, it is impossible to determine 
the sharp value of the tunneling strength, zt for which 
a MI-BG transition can be observed. The variation of 
averaged SF order parameter ip and R for both the odd 
and even MI lobes for two different disorder strengths are 
shown in Fig. [2] at A/C/o=0.3 (dotted (red) lines) and 
A/C/o= 0.5 (dashed (green) lines). 

For the odd MI lobes, we found that the MI region 
starts to shrink due to appearance of the BG phase 
[see Figs. 2(a) and 2(b)] and the region spanned by the 
MI phase gradually decreases at the larger value of the 
disorder strength. When the disorder strength exceeds the 
critical value of disorder for the respective MI lobes, the 
BG phase completely destroys the MI phase, resulting in 
only the existence of the BG and SF phases as pointed 
out earlier. Thus at A/[/o=0.5, the BG region extends 
till zt/Uo ~ 0.072 as perceived from the vanishing of 
the average SF order, ip below 0.1 and finite average 
compressibility, R. 

We also study the behavior of even MI lobes in the 
antiferromagnetic case at different disorder strengths. In 
the pure case, we choose p/Uo=lA and we found that 
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FIG. 2. (Color online) The variation of and at U 2 /U 0 = +0.1 (antiferromagnetic) with disorder strengths A/Uo — 0.3 (dotted 
(red) lines) and 0.5 (dashed (green) lines) respectively and = ij)-, R+ = R- and 'ipo = ^0 ~ 0. The pure case {A/Uo ~ 0) is 
included for comparison (solid (black) lines). At fi/Uo ~ 0.4, ip and R show continuous variations indicating a second order phase 
transition for MI lobes with p = 1 [(a) and (b)]. While at p/Uo = 1.4, ip and R show continuous variation indicating a second order 
transition for MI lobes with p — 2 at A/Uo=0.5. [(c) and (d)]. 


both ip and R show jumps from zero and a constant value 
respectively to distinct finite values signaling the first 
order phase transition for the even MI lobes [see Figs. 2(c) 
and 2(d) (solid (black) lines)]. In presence of disorder, the 
MI state starts to diminish due to the appearance of the 
BG phase and the MI-SF phase transition continues to be 
first order for A/{7o=0.3 [see Figs. 2(c) and 2(d) (dotted 
(red) lines)]. While at A/C/o=0.5 [see Figs. 2(c) and 2(d) 
(dashed (green)lines)], the MI state still survives but the 
MI-SF transition becomes a second order phase transition. 

We have also investigated the behaviour of the individ¬ 
ual SF order parameter and local density components at 
different disorder strengths at both values of chemical po¬ 
tential, namely p/Uq=QA and 1.4. At A/Uo=0, it is found 


that the SF order parameters behave as, ip+='ip-=ipo=0 
with the occupation densities becoming ~ 0.5 and 

Po — 0 at p/Uo=0A and p+=p_=po=2/3 at p/Uo=lA 
in the MI phase. While ip+=ip- ^ 0, ■i/'o=0 and p+, p_, 
po behave similar to that of p in the SF phase satisfying 
the condition p=J2(j Po- correspondingly ipa and Ra at 
^/t^o=0.3, 0.5 are shown in Fig. [2]. The compressibility, 
Ra is constant in the MI phase and gradually increases in 
the SF phase for even MI lobes. 

The above scenario is equivalent to the MI phase with 
fractional occupation densities for an individual spinor 
components, although the summed over spinor components 
yields an integer. 
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FIG. 3. (Color online) The variation of ipa and Ko- at U2/Uo~-0-2 (ferromagnetic) with disorder strengths A/Uo — 0.3 (dotted (red) 
lines) and 0.5 (dashed (green) lines) respectively and k+=k_. The pure case (A/Uo = 0) is included for comparison (solid 

(black) lines). At ^/Uo = 0.4, tp and R show continuous variation indicating a second order phase transition for MI lobes with p = 1 
[(a) and (b)] 


(ii)Ferromagnetic case: For atoms, the exper¬ 

imental data yield 02 < oq for the scattering lengths 
in the respective channels for which the spin dependent 
interaction U 2 < 0 and the spin-1 ultracold atoms show 
similar behavior as that of scalar or spinless Bose gas 
where the MI-SF phase transition is now a second order 
one. 

In the ferromagnetic regime, all the three components of 
the SF order parameter are non zero in the SF phase and 
the ground state is now functions of ipi+, ipi-, ipio- In the 
pure case, we found that ipi^ = tpi- and ipio = V^ipi+ and 
the total spin operator (S)^ in the superfluid phase now 
only takes the value of unity which confirms that we have 
a ferromagnetic superfluid state [19], [20]. 

The variation of ^p and R are shown in Figs. [3(a) and 
3(b)] for both pure and disordered cases at U 2 /U 0 = —0.2 
and h/Uq = 0.4. When A/Uq = 0, the local density, p = 1 
and '0, R are zero till zt/Uo=0.13 (< ztc) indicating the 
system is in the MI regime. As zt/Uo increases, the system 
smoothly goes toward the SF regime with finite ip and 
R and non integer local densities [see Figs. 3(a) and 3(b) 
(solid (black) lines)]. 

As disorder is introduced, the BG phase intervenes in 
between the MI and SF phases. Thus the width of the MI 
phase at a disorder strength, A/C/o=0.5 [see Figs. 3(a) and 
3(b) (dashed (green) lines)] becomes much lesser as com¬ 
pared to that of A/[/o=0.3 [see Figs. 3(a) and 3(b)(dotted 
(red) lines)]. From Fig. [3], we see that at A/C/o=0.5, the 
region spanned by the MI phase is completely occupied 
by the BG phase, resulting in only the survival of the BG 
and the SF phases owing to similar reasons as that of the 
antiferromagnetic case discussed earlier. 

We also study the behavior of the individual SF order 
parameter and local density components at different disor¬ 
der strengths at p/Uo=OA. At A/Uo=0, ip+='ip-='ipo=0 


and p+=p- ~ 0.25 and po ~ 0.5. The compressibility, R^ 
is constant in the MI phase and corresponding ip^ and R^ 
at A/C/o=0.3, 0.5 are shown in Figs.[3(a)and 3(b)]. 

B. Indicators of MI, BG and SF phases 

In the previous section we have studied the effects of 
onsite disorder on the averaged SF order parameter and 
the compressibility. Since we have averaged the SF order 
parameter and the compressibility, it becomes very difficult 
to locate the precise parameter values for a MI-BG and 
BG-SF transition, owing to the fact that the ipi and the 
Pi are explicitly site dependent. Thus at the transition 
point there are some sites with zero SF order parameter 
and integer occupation densities, while other sites may 
have zero SF order parameter and non integer occupation 
densities. 

In order to circumvent this difficulty, we define a 
measurable quantity, y which is defined as. 

Sites with ipi ^ 0 and pi ^ integer 
Total numbers of sites 

In a similar manner, we can also define the Xa- involving 
tpirj and pia- which will aid in the study of the individual 
spinor components. 

The main reason to define x is to distinguish between 
the MI, BG and SF phases where depending upon the 
value of X, we will be able to easily characterize the three 
different types of phases. 

At this point, one can ask, what will be the value of ipi 
and Pi which can define the MI or the SF phase. Since 
we see from the behaviour of ip and R in the pure case 
that the system will remain in the MI phase till zt < ztc 
having vanishingly small tpic- (values below 0(10“^)) and 
the occupation density pi = N ± 6 where N = 0,1, 2,3... 
and i5 is of 0(10“^). These values are set for numerical 
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FIG. 4. The variation of x<t in the antiferromagnetic case with U 2 /Uo=-\-G.l for both the odd and even MI lobes corresponding to 
the pure case A/[/o=0.0 (dotted (black) lines) and for disordered values, A/C/o=0.3 (dashed (red) lines) and 0.5 (solid (green) lines) 
respectively. x<y for odd MI lobe at /r/I/o=0.4 is shown in (a) while for the even MI lobes at ^/Uo=lA in (b). 


convergence of the parameters corresponding to the MI 
phase. Again as earlier, we include discussion of the 
antiferromagnetic and ferromagnetic regimes one by one. 

(i) Antiferromagneic case: For antiferromagnetic 

interaction with t/2/t7o=O.I, we study the the variation of 
X with tunneling strength, zt/Uo in both pure and disor¬ 
dered cases at two different values of fi/Uo corresponding 
to the odd and even MI lobes respectively as shown in 
Fig. [4]. 

In the disorder free case for /i/17o=0.4, we see that 
X makes direct transition from 0 to 1 emphasizing an 
unhindered MI-SF transition. When zt < ztc, the system 
is in MI phase, that is all sites have ipi=0 and pi=integer 
and thus x takes the value 0. While for zt > zt^ the 
system becomes a SF phase with 0 and pi ^ integer 
and thereby x assumes a value 1 [see Fig. 4(a) (dotted 
(black) lines)]. As soon as disorder is included, the direct 
transition of x from 0 to 1 is prohibited indicating a 
gradual increase due to the presence of the BG phase. 

In presence of disorder, that is at A/[/o=0.3 [see Fig. 4(a) 
(dashed (red) lines)], x remains zero till zt/Uo = 0.133, 
beyond which x starts to increase indicating that some 
sites with pi=integer in the MI regime start evolving and 
the system move towards the BG regime with pi ^ integer. 
Further, the gradual increase of the slope of x refers the 
intrusion of the BG phase into a territory which used to 
be the MI phase without disorder. As disorder strength 
increases, that is at A/[/o=0.5 [see Fig. 4(a) (solid (green) 
lines)], we see that x only takes finite values and ultimately 
go unity as a function of zt/Uo, indicating that the region 
spanned by the BG phase continues to increase and results 
in vanishing of the MI phase. Thus at A/[/o=0.5, the 
system only consists of the BG and SF phases as expected 
since the disorder strength A/C/q exceeds the energy gap 


for the first odd MI lobe (p=l). 

It is now quite interesting to study the behavior of Xo- 
incorporating the effects of disorder. In order to study 
Xcr, we use ipia- and pia- and remember that they satisfy 
the same condition as the total SF order parameter, tpi 
and the local density, pi for rendering the MI and the SF 
phases respectively. 

At p/[/o=0.4, the system is in the MI regime since the 
SF order parameter i/i (=\/V^) is small and the local 
densities are as, pi=l, pi+=pi- ~ 0.5, pio — 0(10“^). 
Hence following the definition of Xi we compute x± by 
suitably modifying the numerator of the Eq.(16) with 
tpn- ^ 0 and pi± ^ 0.5 respectively. In Fig. [4(a)], we 
have plotted x± with zt/Uo- We find that x+=X- where 
they show same behaviour as that of x in the clean and 
disordered cases. 

In the antiferromagnetic case, i/io is always zero and pio 
behaves as p± corresponding to the MI and SF regimes 
respectively. So for determining xoj we set the condition 
i/io=t) and pio > 0.001 in the numerator of Eq.(16). 
Its variation with respect to zt/Uo is also shown in the 
Eig.j4(a)]. 

Similar to the discussion above, at /i/C/o=1.4, in the 
MI regime, the local densities yield an even MI lobe with 
Pi=2, pi+= pi-=pio=2/‘i ± 5. So we compute x± by 
replacing the numerator of the Eq.(16) with t/ii± ^ 0 and 
Pi± ^ 0.667 ± (5. For computing xoi we set the condition 
'0io=O and pio ^ 0.667 ± <5 respectively. The variation of 
X and Xa with different disorder strengths are shown in 
Fig. [4(b)]. 

From the behaviour of x, we see that at A/17o=0.3 [see 
Fig. 4(b) (dashed (red) lines)], although there is little value 
in between 0 and 1, but the region occupied by the MI 
phase shrinks as compared with the pure case [see Fig. 4(b) 
(dotted (black) lines)], owing to the emergence of the BG 















FIG. 5. The variation of in fh® ferromagnetic case with 
U 2 /Uo=- 0-2 corresponding to the pure case A/I7o=0.0 (dotted 
(black) lines) and for disordered values, A/Uo=0.3 (dashed (red) 
lines) and 0.5 (solid (green) lines) respectively. 


phase. For A/?7o=0.5 [see Fig. 4(b) (solid (green) lines)], 
the BG phase appears as x takes finite values in between 
0 and 1, while the MI phase still exists since the critical 
disorder strength, Ac/Uq for the p=2 MI lobe is a much 
higher as discussed in the previous section. Components of 
XtT show similar kind of behavior as that of x in presence 
of disorder. 

(ii)Ferromagnetic case: For ferromagnetic interac¬ 
tions with U2/Uo=-0.2 at /i/17o=0.4, the variation of Xcr 
with tunneling strength zt/Uo is shown in Fig. [5]. In 
the ferromagnetic case, since the spin-1 ultracold atoms 
shows similar phase diagram as that of a scalar Bose gas 
(spinless), so there is no distinction between the odd and 
the even MI lobes, and thus showing results involving 
either will suffice. 

At p,/Uq=0A, the MI lobes has local density pi=l and 
X is zero till zt/Uo=0.133. Beyond this value, the system 
goes to the SF phase with x=l in the pure case [see Fig. 5 
(dotted (black) lines)]. In presence of disorder, that is at 
A/C/o=0.3 [see Fig. 5 (dashed (red) lines)], x remains zero 
till zt/C/o=0.085, beyond which y starts to increase and 
at A/[/o=0.5 [see Fig. 5 (solid (green) lines)] y only takes 
the value between 0 and 1 indicating that the system only 
consists of the BG and SF phases. In the latter case, 
disorder completely annihilates the MI phase because 
the disorder strength crosses the critical disorder value, 
AjUo=0A. 

At p/Uo=0A, in the MI regime, the SF order parameter, 
{=y/Wa) is small and the local densities are given 
by pi=\, pi+=pi-=Q.2b ± (5, pio=0.5 ± 5. We compute 
the y±,o by modifying the Eq.(16) with V’i±,o ^ 0 and 
Pi±,o respectively, since in the SF regime, all the three 
SF spinor components are finite. In Fig. [5] we have 
plotted y±,o as a function of zI/Uq which show same be¬ 
haviour as that of y in both the clean and disordered cases. 


C. Percolation analysis and finite size scaling 

At this point we will try to figure out one of our earlier 
questions that was posed in this regard, whether is it 
possible to locate the transition point for the MI-BG 
transition in the presence of disorder with confidence. 
Further what will be the extent of the BG phase, that is, 
whether the whole region is a BG phase or there is certain 
amount of SF phase still left in presence of disorder. To 
answer the questions we have to really look at the concept 
of the appearance of percolating cluster of sites with non 
integer occupation densities and finite SF order parameter. 
In order to understand what do we mean by a percolating 
cluster, we have shown the typical real space plots of 'i/'i 
and Pi for three different values of zt/Uo corresponding 
to the three different phases of the system in Fig. [6]. 
The pair of figures shown in each column depicts MI [see 
Fig. 6(a)], BG [see Fig. 6(b)] and SF [see Fig. 6(c)] phases. 
The top and bottom plots in each of them illustrate the 
values of SF order parameter, x/i and the local density 
Pi respectively in a square lattice of size L x L=50 x 50 
for a single realization of the disorder. The parameter 
corresponding to the plots are A/[/o=0.5 and p/Uo=QA in 
the antiferromagnetic case (t/2/t/o=-l-0.1). The light green 
circles are the indicators for the zero SF order parameter 
and the integer densities, while the deep red circles denote 
non zero SF order parameter and non integer densities. 
The ipi and pi for the MI phase at zt/C/o=O.OOI and for 
the BG phase at zt/t/o=0.06 and for the SF phase at 
zt/t/o=0.3. 

In the clean state, the system splits into two, the MI 
islands with all sites having zero ^pi and integer densities 
Pi as shown by light dots in Fig. [6(a)] till zt/Uo < ztc/Uo 
and the SF islands with finite i/i and non integer densities 
Pi at zt/Uo > ztc/Uo separated by a sharp boundary. As 
soon as disorder is introduced, the BG phase tries to mix 
two islands by removing the boundary between them. So 
the sites with zero i/i and integer densities attempt to 
evolve towards zero tpi and non integer densities as shown 
in Fig. [6(b)] (BG phase), which finally percolates for the 
first time towards a state with non zero i/i and non integer 
densities Fig. [6(c)], the latter being the SF phase. 

To extract different phases of the system, we shall 
concentrate on the appearance of the SF clusters which 
percolate throughout the entire lattice for the first time. 
In the MI states, all sites have integer particle densities 
and hence there will be no SF cluster, while in the BG 
phase, some of the sites have non integer densities and 
thus the SF clusters will be trapped by other clusters with 
sites having integer densities. Finally with the increase 
in the tunneling strength, the SF clusters will start to 
percolate throughout the lattice resulting in a SF phase 
for the system. Thus the BG region can be identified in 
between the MI and the SF phases where SF clusters exist, 
however do not percolate across the lattice. 

When a SF cluster percolates through an infinite lattice 
it is known as an infinite cluster, while when it percolate 
through a finite lattice it is known as a percolating or 
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FIG. 6. (color online) The real space plots of SF order parameter, tpi (first row) and the occupation densities, pi (second row) for 
lattice size L x 71=50 x 50 at A/(7o=0.5 with U 2 /Uo=+Q.l and p/Uo=QA corresponding to three different phases that is the MI, BG 
and SF phases of the system respectively. The light green circles (light shades) represent the zero SF order parameter and integer 
densities and deep red circles (deep shades) represent the non zero SF order parameter and non integer densities. The parameter 
value is at zt/C/o=0.001 in the MI phase (a) and zt/[/o=0.06 in the BG phase (b) and zt/Uo=O.S for the SF phase (c) respectively. 


spanning cluster. Since we are dealing with a finite size of 
the system for our numerical work, we shall aim to find a 
spanning cluster using HK algorithm [30]. 

Whether a SF cluster will percolate can be well under¬ 
stand from the following quantity, Xoo which is defined as 

[40], 

Sites in a spanning cluster 
Total number of occupied sites 

Here we briefly outline the idea used in the HK algorithm 
for the sake of completeness. It is based on the special ap¬ 
plication of the union-find algorithm where it aims to assign 
a label to each cluster. For that, we represent our 2D lat¬ 
tice as a matrix and label all occupied sites initially with 
-1 and unoccupied sites with 0. Our cluster index starts 
with 1. Corresponding to each occupied site, every time 
we check the neighbours at the top and left corner of the 
current site. If both sites are empty, we label a new cluster 
number that has not used so far. Else, if the site has one 
occupied neighbour, then we assign same cluster number to 
the current site. If both neighbouring sites are occupied, 


then we set the smallest number cluster label of the occu¬ 
pied neighbours to use as the label for the current site. To 
link two clusters, we create a union between both labels 
and set the site as the lowest of the two labels. When we 
burn the lattice for a second time, we collect the unions 
and update the lattice. 

In Fig. [7(a)] we study Xoo as a function of x at a disorder 
strength A/C/o=0.5 and /i/C/o=0.4 in the antiferromagnetic 
case. It is found that Xoo is zero till x=Xc=0.581 for a lat¬ 
tice size L X L =128 x 128 which is close to the critical 
threshold value for the occupation probability Pc=0.592 for 
random site percolation problem in an infinite 2D square 
lattice [40] and shows finite value above X > Xc- These 
similarity is only coincidental as the random site percola¬ 
tion model is a classical problem where sites are randomly 
occupied, while we have a system of quantum particles and 
the driving parameters are systematically varied at a given 
value of random disorder. 

Another useful quantity in this regards, is the mean clus- 
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FIG. 7. (Color online) The percolation probability, Xoo and the mean cluster size, Mca corresponding to a lattice size L x L=128 x 128 
at A/Uo=0.5 for the antiferromagnetic case with U2/Uo~+0.1 and fj,/Uo=0.4: are shown in (a) and (b). The percolation probability, 
Pperc corresponding to different lattice size L with system parameter, r]=zt/Uo and finite size scaling are shown in (c) and (d). The 
critical tunneling strength rjc where Ppsrc for different L values intersect at ric=ztc/UQ=0.073. 


ter size, Mcs defined as [40]- 


Mcsix) 


ETp^^pix) 

ETp^pix) 


(17) 


where pSp is the number of occupied sites belonging to a 
p-th cluster and the spanning clusters are excluded from 
the sum. Since the spanning clusters are excluded from the 
sum, the mean cluster size continues to increase with x till 
the appearance of a spanning cluster for the first time and 
starts to decay immediately after that. The variation of 
Mcs with X at a disorder strength A/C/o=0.5 for p/Uo=0.4: 
[see Fig. 7(b)] shows that it reaches its peak value just be¬ 
low the Xc and falls off after Xc which confirms that the 
system is in the BG phase till Xc, and beyond that it goes 


toward the SF phase. 

To deal with the finite size effects and also to determine 
the percolation transition, we resort to the finite size scal¬ 
ing. For a given system size L, we shall now study the per¬ 
colation probability Pperc which is the probability of having 
a percolating cluster with the system parameter r]=zt/Uo- 
Fig. [7(c)] shows the variation of Pperc with different lat¬ 
tice sizes L at a particular disorder strength, A/Uo=0.5 
and p/Uo=QA in the antiferromagnetic case. We averaged 
our results over 1000 (L=10) and 100 (T=30, 50) different 
realizations of disorder. The Pperc is assumed to follow a 
scaling law near the critical tunneling strength, rjc [40] and 
is described by, 


Pperc{L, r]) = p{L - {rj - pc)) 


(18) 
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FIG. 8. (Color online) Phase diagram of spinor ultracold atoms in t — p plane at U 2 /Uo=+Q.l (antiferromagnetic). In the pure 
case (A/[/o=0.0), the odd Ml lobes start to shrink while the even Ml lobes expanded and the colour variation corresponds to the 
magnitude of the SF order parameter ^|) shown in the legend of (a). Phase diagram based on the information obtain from x for the 
disordered cases, namely A/Uo =0.3 is shown in (b) and A/[/o=0.5 in (c). 


where Pperc for different L values intersect approximately 
at one critical point given by r]c=ztc/Uo and p is the scaling 
function which approaches zero in (rj — rjc) << 0 and unity 
for [rj — rjc) >> 0, is the critical exponent which is equal 
to 1.33 for conventional random site percolation problem 
in two dimension. The finite size scaling plot, where all 
data corresponding to different system sizes collapse on a 
single curve, is shown in Fig. [7(d)] where the corresponding 
critical tunneling strength is given by r]c=ztc/Uo=0.073. 
For the sake of brevity we have not included the discussion 
corresponding to the ferromagnetic case. 

We are now in a position to distinguish the three 
different types of phases and able to ascertain the critical 
tunneling strength, r]c=zt/Uo corresponding to the Ml-BG 
and the BG-SF phase transitions. As we have already 
pointed out earlier that the Ml regime corresponds to 
X=0, while the BG regime corresponds to 0 < x < Xc 
till Xoo=0 and the SF regime is given by, Xc < X ^ 1 
when Xoo 7^0. For that purpose we shall tune zt/Uo in a 
controlled manner to closely monitor the value of x and 
Xoo for a given lattice size namely, L=30. We find that 
in the antiferromagnetic case (at ^/Uo=0A), x is zero till 
zt/Uo=0.122 at A/t/o=0.3, while it is non zero at higher 
disorder concentration A/Uq=0.5. Thus at h/Uq =0.4, the 
onset for the Ml-BG transition occurs at zt/UQ=0.122 and 
the BG phase extends upto zt/Uo=0.133, where Xc=0.524 
with Xoo=0. Hence the BG-SF transition takes place at 
zt/Uo=0.l33 corresponding to A/[/o=0.3. For A/[/o=0.5, 
since there is no Ml phase, so the BG-SF phase transition 
occur at zt/17o=0.073 with Xc=0.549. 

D. Phase diagram 

In this section we will now present the phase diagram of 
spinor ultracold atoms int — fi plane for both the pure and 
disorder cases based upon the information obtained from 
ipi, X and Xoo- In the pure case, the phase diagram can 
be easily obtained from the information of the SF order 
parameter, ijji, while in the disordered case, we determine 
our phase diagram depending upon the range of x as dis¬ 
cussed earlier corresponding to both the antiferromagnetic 


and ferromagnetic cases. 

(i) Antiferromagnetic case: In the disorder free case, 
the phase diagram of spinor ultracold atoms in presence 
of antiferromagnetic interaction with U 2 /Uo=+OA is 
shown in Fig. [8(a)]. It shows that the Ml lobes with even 
occupation densities are somewhat enhanced as compared 
to the corresponding lobes with odd occupation densities 
[19, 38]. The Ml-SF phase transition for the even Ml lobes 
is first order in nature, while it is second order for the 
odd MI lobes, as obtained from the variation of the SF 
order parameter and compressibility which was outlined 
in the preceding section. The MI phase with even number 
of bosons per site where each pair of bosons favour the 
formation of singlet parts and likely are to the more stable 
and keep themselves isolated from nearest neighbour 
tunneling. Whereas the MI phase, with odd number of 
bosons per site, although each pair of bosons can make 
singlet states, the remaining one is free to tunnel to nearest 
neighbour sites rendering the odd MI lobes unstable and 
smaller in width as compared to the even MI lobes. 

In presence of disorder, we now try to obtain our phase 
diagram entirely based on the information obtained from 
X. In the previous section, since we distinguished between 
the three different kind of phases depending upon x and 
Xoo- Thus for obtaining the phase diagram, we set x=0 
as the indicator for the MI phase, 0 < x < Xc (=0.524) 
for the BG phase and 0.524 < x A I for the SF phase 
at A/[/o=0.3. Similarly for A/t/o=0.5, where Xc is equal 
to 0.549. The phase diagrams for both the disorder 
strengths are shown in Figs. [8(b) and 8(c)]. At A/?7o=0.3, 
the intervening BG region for the odd MI lobes is more 
noticeable as compared to the even MI lobes [see Fig. 8(b)]. 
At A/t/o=0.5, the BG phase completely destroys the first 
odd MI lobe, but the even MI lobe still exists as expected 
from the earlier discussion [see Fig. 8(c)]. 

(ii) Ferromagnetic case: In Fig. [9(a)], we plot our 
numerically obtained phase diagram corresponding to 
the pure case with a ferromagnetic interaction namely. 
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FIG. 9. (Color online) Phase diagram of spinor ultracold atoms in t — /r plane at 172/[fo=-0.2 (ferromagnetic). In the pure case 
(A/t/o=0.0), the phase diagram is similar to that of scalar Bose gas (a). Phase diagram based on the information obtain from x 
for the disordered cases, namely A/Uo =0.3 is shown in (b) and A/t/o=0.5 in (c). 





A/Uo=0 

A/Uo 

=0.3 

A/Uo 

lO 

d 

II 




MI-SF 

MI-BG 

BG-SF 

MI-BG 

BG-SF 

Interactions 

MI lobes 

p/Uo 

zt/Uo 

zt/Uo 

zt/Uo 

zt/Uo 

zt/Uo 

AF ([/2/Uo=0.1) 

first odd 

0.4 

0.16 

0.122 

0.133 

No 

0.073 

first even 

1.4 

0.26 

0.210 

0.225 

0.123 

0.174 

F {U2/Uo=-0.2) 

first 

0.4 

0.133 

0.085 

0.103 

No 

0.047 


TABLE I. The transition points for the MI-BG and BG-SF phases obtained from x s^nd Xoo for both antiferromagnetic (AF) and 
ferromagnetic (F) interactions are presented in this table. The corresponding parameter values are mentioned in the table. 


U2/Uo=-Q.2. The phases bear similar characteristic fea¬ 
tures as that of scalar or spinless Bose gas. All the three 
SF spinor components are non zero in the SF phase which 
smoothly vary across the MI to SF phases and hence show 
a second order phase transition. In the disordered case, as 
earlier, we determine the phase diagram by invoking the 
critical value of Xc=0.5034 at A/Uo=0.3 and Xc=0.557 
at A/?7o=0.5 corresponding to the onset of the SF phase. 
The phase diagram corresponding to both the disorder 
values are shown in Figs. [9(b) and 9(c)]. For A/[/o=0.5, 
the MI lobes disappear since the disorder strength is higher 
than the critical disorder strength and the system is left 
only with the BG and the SF phases. 


IV. CONCLUSIONS 

We have studied the effect of onsite disorder in a two 
dimensional SBHM for both the antiferromagnetic and fer¬ 
romagnetic interactions. The appearance of the BG phase 
is observed via the average SF order parameter ip and com¬ 
pressibility K. The ipcr are zero while are zero or constant 
in MI phase and gradually increases with disorder yielding 
signature for a BG phase. In the antiferromagnetic case, 
the MI- SF phase transition for odd MI lobes is always sec¬ 
ond order while for even MI lobes, it shows second order 
phase transition with increasing disorder strength. We also 
find that the MI phase exists till the disorder strength is 


below the critical value, while it vanishes above the critical 
value. 

Further the three different types of phases, namely the 
MI, BG and SF phases are identified based on the concept 
of the SF percolating clusters. We define the MI phase 
where no SF cluster exists and the SF phase with at least 
one SF percolating cluster exists. The BG region corre¬ 
sponds in between the MI and SF phase that is the SF 
cluster exists but does not percolate. We locate the tran¬ 
sition point for the MI-BG and the BG-SF phase by cal¬ 
culating the probability of having a percolating cluster Xcx 
(Pperc)- A percolation analysis of SF percolating cluster is 
studied using HK algorithm with different system sizes L 
as a function of tunneling strength, r]=zt/UQ which obey 
finite size scaling law. 

Finally, we obtain our phase diagram based on x for the 
antiferromagnetic and ferromagnetic cases. For that pur¬ 
pose, we determine the critical value of x where a SF clus¬ 
ters percolate first time throughout the lattice and set the 
respective ranges for x for three different phases. Thus a 
reliable enumeration of different phases can be obtained at 
desired parameter values of the SBHM Hamiltonian. Since 
these parameters can be precisely controlled in experiments 
using the properties of the laser fields which superimpose to 
form optical lattices, our study may have significant impact 
on the experimental scenario to determine these phases [41]. 

The transition points for the MI-BG and BG-SF phases 
obtained from x Xoo are presented in the table I for 
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both antiferromagnetic (AF) and ferromagnetic (F) inter¬ 
actions in pure and disordered cases. 

Now we shall include some comments concerning the 
drawbacks in the mean field theory. Since our work is 
mainly a focus on the effects of disorder in SBHM, thus 
the single site mean field theory always is not able to han¬ 
dle the site inhomogeneity of the system properly and it 
works well in higher dimension d > 2. Although we take 
into the effect the site inhomogeneity by defining but 
ignoring quadratic fluctuations introduces a tiny error in 
the calculations. For example, the MI phase will vanish 


above the critical disorder strength, our results are in good 
agreement against this claim within a 3% error. 
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